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Abstract 



A quantum simulation of a travelling salesman is described. A vector space for a graph is defined 
together with a sequence of operators which transform a special initial state into a superposition states 
representing Hamiltonian tours. The quantum amplitude for any tour is a function of the classical cost 
of travelling along the edges in that tour. Tours with the largest quantum amplitude may be different 
than those with the smallest classically-computed cost. 

1 Introduction 

In the problem of the travelling salesman we are given a graph with a fixed number of vertices. Every vertex 
is connected to every other vertex by one edge. A vertex is not connected to itself. Every edge has a cost 
which is a positive real number. A Hamiltonian tour of such a graph is a path which begins at a special 
vertex, chosen arbitrarily but then fixed. Starting from this vertex, a sequence of steps is taken in each of 
which a previously unvisited vertex is visited by travelling along the edge connecting the two. The cost of 
travelling along that edge is accumulated. The classical cost of a tour is the sum of the costs associated with 
the edges traversed during the tour. The tour ends when all vertices have been visited exactly once at which 
point one more step is taken to return to the starting vertex. For a graph with N vertices, there are (N — 1)! 
tours that begin and end at the starting vertex. The problem is to find the tour with the minimum cost. 
This is the simplest version of the travelling salesman problem pQ. 

In section 2 I will describe a quantum simulation of this problem. The simulation occurs in a vector space 
for a graph. This vector space is sufficient to describe Hamilitonian tours which are special basis states. The 
simulation begins by placing the graph in a basis state corresponding to a choice of the starting vertex. A 
sequence of operators will transform this initial state to a final state which is a superposition of only the 
Hamiltonian tours. In this final state all tours have the same quantum amplitude. There are three goals for 
the description in section 2. The first is to describe the vector space for a graph. The second is to show how 
this vector space is capable of describing a Hamiltonian tour. The last goal is to construct operators which 
perform the transformations which may be interpreted as carrying out all Hamiltonian tours. 

The simulation is improved in section 3 by taking into account the cost of travelling along the edges. The 
operators described in section 2 are modified so that the quantum amplitude of a path is built up as the tour 
proceeds. The amplitude is now not a constant, but a function of the cost along the edges. Hamiltonian 
tours with smaller classical costs have larger quantum amplitudes. 

The dynamical system described is not a quantum algorithm that can execute on a universal quantum 
Turing machine described by Deutsch [2]; it is a simulation of a special system by a quantum mechanical one. 
Such sytems [3 |4l [5] have long been in existence. The simulation requires a quantum mechanical computer 
capable of executing it. The model of a computer most suitable has been described by Feynman [4] based on 
work of Benioff [3] and othertQ. 

Feynman describes a universal quantum computer which is a one dimensional lattice of spin-half systems 
with nearest neighbour interactions coupled to a system of logic gates. The state of the complete system 
consists of the state of the spin lattice and the state of the registers on which the gates act. The Hamiltonian 
for the system defines the dynamics of the lattice of spins as well as that of the gates. Nothing in Feynman's 
description requires the system to be a set of logic gates; they can be more general [5], In the simulation 
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of the travelling salesman presented here, the state space of the registers are replaced by the state space 
of a graph while the operators which act on the registers are replaced by the operators for carrying out a 
Hamiltonian tour. The system coupled to the lattice of spins is not a set of logic gates but one which operates 
on a graph. 

2 Simulating Hamiltonian tours 

Let G be a fully connected, symmetric graph with N vertices. These two properties mean that every vertex 
is connected to every other vertex by an edge and that the same edge is used when travelling from one vertex 
to another in either direction. A vertex is not connected to itself. The vertices are labelled 1, 2, . . . , N. V% is 
the name for the vertex labelled i. Nothing that follows depends on the choice of labelling. 

Associated with G is the quantum system Q which is the system that will simulate the travelling salesman. 
In the following, unless explicitly noted, the graph in question will be the quantum system Q and not G because 
it is Q that is subject to dynamical evolution. 

I will define a vector space for Q and assign meaning to these states. These states will be sufficient to 
define a visit of a vertex and a Hamiltonian tour. Q is prepared to be in a special initial state corresponding 
to starting at the chosen initial vertex. A sequence of N operators applied to this state will transform it to a 
final state which is a superposition of all Hamiltonian tours and only the Hamiltonian tours. Let the initial 
and final states be |<?(0)) and \g(N)) respectively and let S(l), S(2), . . . , S(N), be the operators. The goal is 
to show that 

\g(N))=S(N)S(N-l)---S(l)\g(0)), 

and that \g(N)) is a superposition of states which may be interpreted as Hamiltonian tours. Moreover \g(N)) 
contains only such states and no others. 

Q is driven from |g(0)) to \g{N)) by being coupled to an external system. In other words, the step 
operators S(l), S(2), . . . , S(N) are components in a larger Hamiltonian for the dynamical system consisting 
a one-dimensional lattice of spins and Q [3] . 

The vector space for Q is defined to be the direct product of the space of each of its vertices. The space 
at each vertex is the direct product of (N + 1) qubit spaces. At each Vi, the basis vectors are 

\vi;an,ai2, ■ ■ ■ ,<Hn) ■ (1) 

The values of U{, and ay can be or 1. The basis vectors span a space of dimension 2^ JV+1 - ) giving Q a 
dimension of 2 N ( N+1 \ These states can be transformed by the usual raising and lowering operators. So 

«i|0)=0, «i|l) = |0), t>f|0) = |l), «?|1) = 0, i = l,2,...,N. 
Similarly for the states labelled by a^. 

a 4 ,|0)=0, a ij \l) = \0), a* j \0) = \l), a*.|l>=0, i, j = 1, 2, . . . , N. 

These operators obey the well known commutation rules. 

ViV* +v*Vi = 1, 

ViVj - v*Vi = 0, (i ^ j) 

& ij°-*M - Oij°>kl = 0, ^ (k, I). 

Where k, I = 1, 2, . . . , N. The v and v* operators commute with the a and a*. The basis vectors of Q are 

N 

Y[\vi;aa,a i2 ,...,a iN ) . (2) 
»=i 
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The meaning of Vi, and follows, vi = 1 means Vi is being visited while Vi = means V, is not 

being visited. ay = 1 means there is a directed edge from Vi to V, while ay =0 means there is no such edge. 
A state with an = 1 means there is a self-edge from Vi to itself. For any i the set of values ay represents the 
edge-state of V^. 

Let us first see how these states can describe a Hamiltonian tour. Consider a graph G^ 3 \ with three 
vertices V±, V 2 and V3. The basis vectors of its quantum equivalent are 

flu, 012013) |«2 ; 021 , 022) «23> I v 3 ; a 3i , a 32 , a 33 ) 

for all values of the Vi and ay. Begin by placing its C/'- 3 -' in state in which one vertex V%, is visited and all 
three vertices have only self-edges. According to the convention, in this state 



V\ =1,V 2 = V3 = 0, 

where i, j = 1,2, 3. This is the state 

| 5 ^(0)) = |1; 1,0,0)|0; 0,1,0)|0; 0,0,1). 

A shorthand for this state is 

\g^(0)) = 1 1; au = l) |0; 022^1) |0; a 33 = l ) . 

The underline indicates that all ay not appearing in the ket have a value of 0. In this state, only V\ = 1 so 
V\ is the vertex being visited. All ay = 0, (i 7^ j) and an = 1 so all vertices have only self-edges. The graph 
representing this state is shown in Figure Q] 
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Figure 1: A graph representing the state \g^(0)}. The • marks the vertex being visited, here Vi. Every 
vertex is connected only to itself with a self edge. 

Let us take the first step by acting on |</ 3 '(0)) with the operator uiana^W^- 

viaiia* 12 V2 1 1; on = l) |0; a 2 2 = l ) |0; o33 = l ) = |0; ai2 = l ) |l; a22 = l ) |0; o33 = l ) ■ 

In this new state, v± = and v 2 = 1 so by convention this means V 2 is being visited. The other change is 
that an = while ai 2 = 1 which means there is a directed edge from V\ to V 2 - The effect of the operator 
v\a\\a* 12 v 2 is to make V% the vertex being visited and replace the self edge at V\ with a directed edge from 
V\ to V%. The state of V 3 is unchanged. The meaning of this transformation is that we have carried out the 
first step of one tour from V\ to V 2 - Figure [2] shows the graph representing this state. 
In the second step, visit V 3 from V 2 with the operator £2022023^3- 

(ha^a^v^iv^nal^) |5 (3) (0)) = |Q; ai 2 = l ) |0; a 23 = l ) |l; a 33 = ) 

After this second step (Figure [3]) , V3 is visited and we have recorded that the visit of V3 took place along the 
(2,3) edge of V 2 (a 2 3 — 1). No unvisited vertices remain, so we must return to the starting vertex with the 
operator vsa^d^vl. This produces the state 
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Figure 2: A graph representing the |0; an = 1) 1 1; a 2 2 = 1) |0; 033 = 1) produced as a result of the first step from 
Vi to V2. V\ has no self-edge. V2 is is being visited. 




Figure 3: A graph representing the state |0; 012 = 1 ) |0; o23 = l ) [1; 033 = ). V3 is is being visited. 



{v^aasa^v^i^^a^v^iviaiianv^) \g (3) (°)) = |l; oia = l) |0; a 23 = l ) |0; a 3 i = l ) 

which is a basis state corresponding to the Hamiltonian tour V\ —} V2 — >• V3 V\ . It is carried out with the 
three operators v\a\\a\^v\^ £2022023^3 and v%aiza%\v\. 

This three- vertex graph has only one other Hamiltonian tour — V\ — > V3 — > V% — > V\ — which can be 
navigated with the operators 010110131)3, 030330321)2 an d ^022021^1 • ^ we can perform both transformations 
simultaneously the final state of the graph, ignoring normalisation, would be 

Vi->v 2 ->v 3 ->-v 1 v 1 ^v 3 ^v 2 ^v 1 

|l; oi2 = l) jO; 023 = 1 ) 10; 031 = 1) + |l; oi3 = l ) j0; 032 = 1 ) 10; Q2i = l) 

In both terms V\ is the vertex being visited and, since no vertex has a self edge, all vertices have been visited. 
This state is shown in Figured] This example has shown that the basis states in (2) are capable of encoding 



,2, 
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1. "^^ 3 

Figure 4: The final state |l; oi 2 = l) |0; a 2 3 = l ) |0; a 3 i = l ) + |l; ai3 = l ) |0; 032 = 1) |0; o2i = l ). 
a Hamiltonian tour. 

A Hamiltonian tour starting at V\ is a basis state with the following properties (i, j = 1,2,... , N). 

1. Vi=l J » i = Q > (*^l), 

2. o„ = 0, 

3. for each i, one ay = 1; all other dij =0 (i 7^ j) and 

4. the directed graph represented by the basis state is connected. 



In the general case Q is prepared in an initial state in which V\ is visited and all vertices are in a state 
corresponding to having only self edges. 

N 

|ff(0)) = 1 1; an = l) f[ |0; ajj = l ) . (3) 

i'=2 

To take the first step, define the step operator 5(1). 

1 N 

5(1) = , N _ 1 51 A^v.auatjV* (4) 

i,j'=l 

where Ay is the adjacency matrix of G. 



Aij 



if i=j; 

1 ifi^j- 



Let \g(l)) be the state produced by acting on \g(0)) with 5(1). 

^ N N 

1.9(1)) - 5(1)1.9(0)) = -== Y, A^aua^v* |l; o u = l ) J] |0; o fcfc = l ) 

^ i,j=l fe=2 

In the sum over i the only term that survives is the one with i = 1 because v\ = 1 and i>; = 0, (i 7^ 1). So we 
have 

^ AT N 

1.9(1)) = -7==== J] AxjVihualM |l; o X i = l ) TT |Q; q fcfc = l ) 

^-1^ fe = 2 

Of the remaining terms consider the first with j = 1 which is -An'Oidnd^fiJ' . This does not contribute 
because An = 0. Therefore |g(l)) is a superposition of (JV — 1) basis states and does not contain a state 
resulting from visit of V\ to itself. After the first step we have 

1 N 

|.9(1)) = 5(1)|.9(0)) = -j== Y\l^i*). (5) 

This is a convenient new notation to represent incomplete tours. Suppose we want to represent the state 
for the path Vi — >■ V2 —> V3. Here, V3 is being visited from V\ via V2. None of the other vertices have been 
visited and the tour is currently stopped at V3. That is, no steps whatsover have been taken from V3. We 
will write this as 



|1 -*-2-> 3«). 

The number appearing last in the sequence identifies the vertex being visited, in this example V3 . All vertices 
whose labels appear in the sequence are in a state that is consistent with this path. The bullet after the last 
number signifies that the remaining vertices, whose labels do not appear in the sequence, have not taken part 
in the tour. The tour is stopped at the vertex whose number appears last in the sequence. Written out fully, 

N 

|1 2 -> 3 •) d = |0; oi2 = l ) |0; a 23 = l ) |l; a 33 = l ) II \°'> akk = l J ■ 

In each of the (JV — 1) terms in equation a different vertex is visited. |g(l)) is a normalised state by 
construction. However, 5(1) is not a unitary operator: 

5(1)*5(1) ^ t 

But it can be shown that 
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S(iy\g(l))=S(irS(l)\g(0)) 

= b(o)). 

Equation (jSJ) shows that 5(1)* reverses the action of 5(1) when Q is initially in the state \g(0))- Although 
5(1) and 5(1)* are not inverses of each other, they behave as inverses when they act on \g(0)). 

For the second step of the tour we define the step operator 5(2) which will operate on \g(l)} to produce 
1.9(2)). 

| 3 (2)) =5(2)|. 9 (1)), 

where 

1 N 

= /N-9 ^ A -J r '"""'.i''j.i"jJ r 'r 

5(2) is identical to 5(1) except that it contains a new operator combination &jj<ijj and a normalisation factor 

of (N — 2) not (N — 1) . The role of the additional operator combination is described presently. 5(2) 
will act on the (N — 1) terms of |g(l)) in equation <j5j> to produce many more terms. Let us look at 5(2) 
acting on just the first term of |<?(1)), namely |1 — > 2 •). In the sum over i the only term that survives is the 
one with i = 2 as all others are annihilated by the action of Vi. From N 2 terms we are reduced to N. 

1 N 

5(2) |1 -> 2 .) = -=== J2 A 2j v 2 a 22 a* 2j a* jarj v* |1 -> 2 •) . 

i=i 

The term with j = 2 does not contribute because A 22 =0 so the sum over j produces at most (N — 1) terms. 
The term with j = l corresponds to visiting V\ from V 2 which would certainly not be a legal step since V\ , 
being the starting vertex, has already been visited. But the operator for this step contains o*jfin which will 
produce because in this state an = 0. So the term with j = l does not survive meaning that V\ is not visited 
from V 2 in the second step. So the sum over j produces just (N — 2) terms. Since there is nothing special 
about V 2 , the same argument holds for all other vertices. Therefore, from each term of |ff(l)) in equation ([5]), 
the action of 5(2) will produce (TV — 2) terms to give 

N N 

\g(2)) =5(2)5(l)|g(0)) = , = Y Y \1 -> i -» j •) . 

\g{2)) is normalised. The operator 5(2)* reverts \g{2)) to \g(l)). 



1 1 N 

S(2TS(2 M 2))= g (j V- 2 )|l-,i. 



1 N 

It is convenient then, to define the adjacency operator Ajj(T) 

4 (T) - -^=A ljau a* 3 a* ja33l T = 1, 2, . . . , (JV-1) (7) 

where i = 1, 2, . . . , N, and from which we can define a sequence of step operators, 5(T), for the N steps of 
the tour. 

f ES =1 MyCT)*i, if T = 1, 2, . . . , (N - 1) 
S(T) = \ ^ (8) 

lEj=l"i a i»*fl' 5 l> if T=N. 
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The operator for the last step, S(N), does not produce new terms. Its effect is to visit V\ for the second 
time. After the action of S(N) every term in the superposition will have v% = 1. Starting from the initial 
state |5(0)) we produce a final state \g(N)) where 

\g(N)) = S(N)S(N-l)---S(l)\g(0)). (9) 

\g(N)) is a superposition of (N— 1)! basis states and is normalised by construction. Each of these corresponds 
to a distinct Hamiltonian tour. All terms in the superposition have the same probability amplitude of 
1/V^1)!. 

As has been noted, S(T) is not a unitary operator. Although the operators used in Feynman's description 
are unitary, the following weaker condition is sufficient. 

S(T)*S(T)\g(T-l)) = \g(T-l)) 
This implies that the computation is meaningful only starting from the special state \g(0)). 



3 The cost of a tour 

In the simulation just described all tours have the same quantum amplitude. If a measurement is made on 
\g(N)), one tour is as likely as another. This can be improved by taking into account the cost of a tour. Let 
Wij be the weight of the edge defined as follows. 

Wij = if i = j; < Wij < 1 if i ^ j. (10) 
The Wij are non-negative real numbers which are normalised. 

N 

5>£ = i, (?; = i,2,...,a). (11) 

The w^ are taken to be symmetric: w^ — Wji. We can always impose on the w^ the N conditions in 
equation (TTTj) because there are N(N — l)/2 w^. In equation ([7]) we defined (1) the operator for the first 
step 

^ij(l) = ^/ jy ^ Aijaua^a jj hjj , i,j = 1, 2, . . . , N. 
We replace this with a new definition 

Ay(l) "'<."./'','/'.;/'.;./• ( 12 ) 

For the first step the operator 5(1) will use this new Aij(l) rather than the one defined in equation ([7|). 
Acting on |g(0)) with the new Ajj(l) gives 

N 



iff(i)>=&i0(o)> = 5> M |i->*' 



i=2 

Because of the condition in equation (TTT1) . |g(l)) is normalised and it is still the case that 

S{\g{l)) = StS x \g{Q)) = | 5 (0)>. 

The amplitude for the step V% Vt, which was previously (N — 1) , is now wu. For the second step, 
suppose we could define an operator that is identical to Aij(l) in equation ([T2l . 

i ii (2)=i ij (l) 

and operate on |g(l)). As before |g(2)) will contain (N — 1)(A — 2) terms because we have changed only 
the amplitudes and not the dynamics. After the second step the path V% — > Vi 2 will have quantum 
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amplitude wu 1 Wi 1 i 2 . But recall that in the second step, starting from Vi x , we do not visit V± or itself. 
That is, only (N— 2) other vertices are visited not (N— 1). In any step, starting from a given vertex, the total 
probablity of all terms produced must be unity. The total probability of all terms produced by travelling 
from Vi ± to (N — 2) other vertices is produced by summing over z 2 — the index of the vertices visited from 

N 

Therefore |<?(2)) = S(2)\g(l)) is not a normalised state. We are missing a probability of wf il . We can get 
around this situation by using a modified operator for the second step 

An (2) = a H a*^a 33 . (13) 

This is exactly the same as in equation (IT21 except that all weights have been rescaled. In terms 

of the rescaled weights, in the second step, the total probability of all terms produced by starting from Vi 1 
is now 



N 2 N 
sp W M 2 = 1 V- 2 

= 1 



as required. The probability of all paths starting from V± is produced by summing over i\ and %2- 



N 2 N 

E W lii \ ~< 2 

i l= l ^ W ^ l2 ^l 
N 

= £<• 

n=i 

This produces, from equation (|11[) . unity proving that |g(2)) is normalised. 

The rescaling of weights works in step 2 because in this step, every vertex, Vi 1 , visited in the first step was 
visited from the fixed vertex V\. At the end of the first step, every path has the same history — all start at V±. 
This is not the case in subsequent steps. In the third step look for example, at the path V\ — > V% — > V3 — > V4. 
In going from V3 to V4, neither V\ nor V2 is visited from V3, the missing probablity is (tUoj + w 2 2 ) which 
means all amplitudes for the step starting from V3 must be scaled by (1 — (tufj + w^))' 1 ^ 2 - Now look at 
the path V\ — > V3 — > V4 — > V§. Here the missing probability in going from V4 to Vq is (w^ + wf 3 ). So for 
this path the scaling factor is (1 — (w^ + w| 3 )) -1 / 2 . For the second and subsequent steps every path has a 
distinct history which cannot be accounted for by the technique used in step 2. 

We need an operator to rescale the weights. Let us look at this one. 

N 

»•< y2"V'.""'r » 1 1 : ' 

This is a real, linear operator which has every basis state of G as an eigenvector. Let us apply W3 to the path 

v 1 ^v 2 ^v 3 . 

w 3 |1 2 -> 3 •> = (it; 2 ! + wj 2 ) |1 ->■ 2 ->• 3 •) . 

So W3 picks out the total probability along all edges from V3 leading vertices already visited along the path 
Vi V2 — >• V3. In general, if along some path, Vi is the vertex being visited and we operate on that state 
with Wi, the eigenvalue will be the sum of the squares of weights from Vi to all vertices which have already 



N 2 

V w 2 Wlli2 - 

2^ W l*i (1 _ ,,.2 ) - 
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been visited along the path. However the operator we really want is one related to u>i . Look at the rescaling 
operator 



which we can write as a power series 

N N 2 

Wi = c + ci ^2 w %m a n + c 2 (X! "'v ".;,■"./,•) + • • ■ ■ 

4=1 j=i 

where the {cq,c\, . . .) are real numbers. This is a power series of commuting operators a^&n. So Wi is a 
linear operator and can be used to renormalize weights along the edges. The problem is that the Wi are 
highly degenerate and have a large number of eigenvalues which are infinite. Any state corresponding to a 
Hamiltonian tour — in which all vertices have been visited — will be a state with an infinite eigenvalue. 
A Hamiltonian tour is a state in which cijj = for every j. In this state, for any fixed i and distinct 
it, i 2 , . . . ,i(jv-i), 



N N 

(1 - ^ w ii°-3j&jj) I 1 ^ *1 ^ > i(N-l) -> 1} = (1 - X] ™!) I 1 ll ~* > HN-l) 

3=1 3=1 

= 0. 



To avoid this we redefine the scaling operator. 



X -1/2 



W i =(l-v 1 v* 1 ^2wla* jj a jj ) . (15) 

3=1 

The infinite eigenvalues are avoided because the action of v\v\ in equation (|15p will produce for any state 
in which v\ = 1. At both the start and end of a tour the second term in equation (|15[) vanishes to produce a 
scaling of unity which is exactly what is required. 
The adjacency operator for the graph is now 



Aij = Wijaua^&^cijjWi. 



The step operators are 



f £U=i iiA » {T ^h if ^ = 1, 2, ... (AT - 1) 
S(T) = \ ■ (16) 

These will produce all Hamiltonian tours as desired. The amplitude for the last step is always unity. This is 
because in the last step we must return to the starting vertex. But in fact it is unity even in step (A — I). 
Scaling increases the weights along the edges which is expected since after each step there are fewer edges 
along which to leave a given vertex. The amplitude for the tour V\ — > — > Vi 2 — > Vi 3 — > ■ ■ ■ — > Vi N — > V\ 
will be 

w ln , Wtlt2 . Wl213 (17) 
^-<i^-«i+<n) 

where each term in the product is the factor produced in each step. Each factor is a ratio that approaches 
unity as the tour progresses. 

To connect with the classical result, let c^- be the classical cost of travelling along the edge Each 
Cij is a positive real number with = c,, and Cu = 0. Let 
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N 

(C0 2 = £ e~ 2c -- (18) 



An obvious choice for is 



These lUy satisfy the condition in equation (|lip and can be treated as quantum amplitudes. Smaller classical 
costs Cjj will produce larger quantum amplitudes for a step along the corresponding edge. 
In terms of the Cjj equation (|17[) is 



But for scaling, which is an effect of the simulation, the quantum amplitude is the exponential of the sum of 
the classical costs of the edges in the tour. 

The scaling factors have two quite different effects. In the classical analysis, the Hamiltonian tour V\ — > 
V2 — > V3 • • • — > Vn — > Vi has the same cost as the tour in the reverse direction. In this quantum simulation 
these two tours will generally have different quantum amplitudes because the scaling factors will be different 
for each. 

In the classical case the optimal paths are those with the smallest total cost. In this simulation the 
optimal paths are those with the largest quantum amplitude. Due to scaling, Hamiltonian tours with the 
largest quantum amplitudes may not be the same as the classical tours with the smallest classical cost. The 
two results are expected to be the same in the majority of cases but this is based on the general notion that 
the classical optimal path will play the leading role even in the quantum domain. Numerical simulations on 
small problems and a few choices of edge-costs seem to support this claim. 



4 Summary 

By doing quantum mechanics on a graph I have described a simulation which can efficiently traverse all 
Hamiltonian tours and only the Hamiltonian tours. It is expected that in the majority of cases, tours with 
the smallest classical cost will have the largest quantum amplitude. There will certainly be cases in which 
the two differ. One can take this to mean that this simulation does not always produce a final superposition 
containing the correct — classically defined — results. On the other hand the final state contains basis states 
which represent the quantum optimal paths. Computing either is equally difficult. 

It will take 0(N(N — l)/2) measurements on the final state to reduce it to a basis state corresponding 
to a Hamiltonian tour. In general the amplitude for any give tour is 0(1/ y/ (N — 1)1 so as expected, the 
simulation will have to be executed many times to find the optimal paths. This simulation will therefore be 
effective in those cases where the optimal tours have much larger quantum amplitudes than the sub-optimal 
tours. 

The scaling operators Wi have additional infinities. States with V\ = Q and all an = are states with an 
infinite eigenvalue. However these states are never reached in this simulation. 

The marking scheme described in section 2 is by no means unique. There is plenty of room for different 
choices of basis states and the step operators. The one described here is natural because of the central role 
of the adjacency matrix. It is also one of the simplest. The marking scheme described in section 2 is an 
efficient way to traverse all Hamiltonian tours without any scaling. 

One of the objectives of this simulation was to make quantum amplitudes depend in some way on the 
classical cost of travelling along an edge. Other simulations are possible in which the classical costs are used 
in different ways. 
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